d/dx (x (x^2+1)^n)